Spin dynamics in the strong spin-orbit coupling regime 
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We study the spin dynamics in a high-mobility two dimensional electron gas (2DEG) with generic 
spin-orbit interactions (SOIs). We derive a set of spin dynamic equations which capture the purely 
exponential to the damped oscillatory spin evolution modes observed in different regimes of SOI 
strength. Hence we provide a full treatment of the D'yakonov-Perel's mechanism by using the 
microscopic linear response theory from the weak to the strong SOI limit. We show that the 
damped oscillatory modes appear when the electron scattering time is larger than half of the spin 
precession time due to the SOI, in agreement with recent observations. We propose a new way to 
measure the scattering time and the relative strength of Rashba and linear Dresselhaus SOIs based 
on these modes and optical grating experiments. We discuss the physical interpretation of each of 
these modes in the context of Rabi oscillation. 



I. INTRODUCTION 



In recent years research in semiconductor based devices 
has incorporated the spin degree of freedom as a new 
state variable in novel electronic devices with potential 
for future applications. The SOI is a key tool to electri- 
cally manipulate the spin and realize such devices. How- 
ever, the SOI is a double-edged sword because it will also 
induce random spin precession through an angle tt so T be- 
tween collisions with impurities, where r is the electron 
life time. This is known as the D'yakonov-Perel's(DP) 
mechanismir.2. and dominates the spin relaxation in the 
technologically important III-V semiconductors^ There- 
fore it is very important to understand fully the DP 
mechanism for the possible application and further devel- 
opment of spintronics devises. Although the study of DP 
mechanism in semiconductors in the presence of SOI was 
initiated long ago, most of the theoretical researcher— 
focuses on the weak spin-orbit coupling (SOC) regime 
where fl so T -C 1. However, as high- mobility 2DEG sys- 
tems are created, it is now not difficult to reach the strong 
SOI regime experimentally where VL so t > 1 at low tem- 
peratures as long as the mobility is approximately larger 
than 1.2 x 10 5 cm 2 /Vsii 2 - The spin evolution in this regime 
is observed to be damped oscillations in the unifor m 11 ' 12 
and nonuniform spin polarized system^ - — which can 
not be described by spin-charge drift-diffusion equations 
derived for the weak SOC regime and lacks a clear theo- 
retical explanation. 

Here, we study the spin dynamics theoretically from 
the weak to strong SOC regime. The method we use is 
linear response theory 5 ' 7 ! 16 We derive a set of spin dy- 
namical equations in the uniform spin polarized 2DEG 
with different SOIs. We show analytically that for 
fl so T > \ the damped oscillations appear. The decay 
rate in this case is proportional to — instead of r as in 
the weak SOC regime. The cubic Dresselhaus term is 
shown to reduce the oscillatory frequency and increase 
the decay rate in the strong SOC regime. The spin dy- 
namics for non-uniform spin polarization with spatial fre- 



quency q in the strong SOC regime is obtained by solving 
the equations numerically. We discuss these dynamics by 
using the analogy with Rabi oscillations between two mo- 
mentum states which are gaped by the SOI. Our results 
match the experimental observation quantitatively. We 
also show how to exploit our analysis to create an accu- 
rate measurement of the strength of Rashba and linear 
Dresselhaus SOIs in a 2DEG, hence allowing a full char- 
acterization of different device samples which will lead 
to a more accurate modeling and predictability of the 
optimal operating physical regimes. 



II. MODEL HAMILTONIAN AND DENSITY 
MATRIX RESPONSE FUNCTION 

Normally in the 2-D semiconductor heterostructures, 
we have three kinds of SOIs, namely the linear 
Rashb a 19 ' 20 term and the linear and cubic Dresselhaus 2 ! 
terms. The Hamiltonian takes the form 



H = 



*1 
2m 



h(k) 



(1) 



where h(k) is the effective magnetic and contains Rashba, 
linear and cubic Dresselhaus terms which are 



ky 7 k x ) j 



h fl (k) = a( 
h Dl (X) = Pi{k y , k x ), 
h D3 (k) =-2/3 3 cos20(-fc !/ ,fc x ), 



(2) 
(3) 
(4) 



where kf is the Fermi wave vector. Here we take 9 as the 
angle between the wave vector k and the [110] direction 
which is the x axis in our coordinates. The above SOIs 
split the spin-degenerate bands and dominate the spin 
dynamics in the 2DEG. The corresponding SOC Hamil- 
tonian and the spin precession frequency fl so takes the 
form: 



(Ai - 2/3 3 cos 26)k x a y + (A 2 + 2/3 3 cos 29)k y a x , 



(5) 



where Ai = a + /3i, A2 = /?i 
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We derive the spin dynamic equations from the density 
matrix response function^ used previously in Refs. [5| 
and0 in the weak spin-orbit coupling regime fl so T <C 1. 
The spin diffusion is dominated by the pole of the spin- 
charge diffusion propagator or " diffuson" : 

v = [i-i}- 1 (6) 

and 

f f cPk A R 
I(tict2,ct 3 <j4, = J ■^p-G (T3Cri (fc,0)G (T2Cri (fc + g, fi), 

(7) 

where Ui is just a number which can be 1 or 2J± Here, 
we will calculate the diffuson matrix exactly and make it 
valid in any SOC regime. It is more convenient to write 
the Eq. [7J in a classical charge-spin space 

I a0 = Tr(aja p ), (8) 

where a, /3 = c, x, y, z.— 



III. UNIFORM SPIN POLARIZATION 

In the case of a uniform spin polarized 2DEG system, 
i.e. q = 0, because the effective magnetic field due to the 
SOI has inversion symmetry in momentum space, only 
the diagonal elements of the diffuson matrix are nonzero, 
which means the spin x, y, z and charge are not coupled 
to each other. Therefore, when considering the uniform 
spin polarization along the z direction, only I zz needs 
to be calculated. First, we neglect the cubic Dressel- 
haus term which is normally much smaller than the linear 
Dresselhaus term. We find the pole of diffusion matrix 
by solving the equation 

i-i» = i- , l - mT =0, 

v/((l - *nr)2 + (ft so r)2)2 _ 7 2(n soT )4 

(9) 

where fl is the frequency of the spin evolution, f2 so = 

2y/a 2 + Plk f , 7 = = k f is the Fermi wave 

vector and I zz is obtained from the exact angular integra- 
tion of Eq. ([7]). The details of calculating I zz are shown 
in the appendix. The eigenmode of the spin dynamical 
evolution takes the form 



(10) 



f 2 - V2^ i - 2(o so t) 2 + y/i - A{n soT y + 4(n so ^A 



Note that 7 < 1 and 

(1 - 4(O so r) 2 + 4(fi so T)V) < (1 - 2n 2 so r 2 ). 

Therefore, as long as 1 - 4(ft so r) 2 + 4(fi so r) 4 7 2 < 0, 
a non-equilibrium spin polarization will exhibit damped 
oscillation with respect to time. 

In the case of a — or /?i =0, the eigenmode takes 
the form 

tilT = 1-1^13402^2. (11) 

When f2 so r > 1/2, the decay rate changes from the ex- 
ponential decay mode to the damped oscillation mode. 
The oscillatory frequency in the clean limit, r — > 00, is 
f2 so . Several experiments^— observe the damped oscil- 
lation mode of spin evolution at low temperature. How- 
ever their analysis did not explain quantitatively when 
this kind of mode appears but just qualitatively argue 
that it appears in the regime where Q so t > 1. Our the- 
ory agrees with the recent experimenljii in which the au- 
thors observe that when the temperature is above 5 K, 
the oscillation will disappear. In their system this cor- 
responds to f2 so T* rj 0.48, which is close to our result 
1/2. Here t* is different to the transport scattering time 
t p obtained from the mobility; this difference is due to 
the Coulomb interaction effect on spin-currents and spin 



dephasing. 8,13 This e-e interactions treatment is beyond 
our paper and will not be discussed in this work. The r 
here is corresponding to r*. When the oscillatory mode 
appears, the damped decay rate is always equal to J- 
when either a — or j3\ — 0. This result matches the re- 
cent experiment^ in which the authors found the decay 
rates for several different 2DEGs always equals j-j^r when 
the damped oscillatory mode appears, in agreement with 
our theoretical result. 

As the linear and cubic Dresselhaus terms always co- 
exist, we have to consider the effect of cubic Dresselhaus 
term on Eq (fTTj) . We do this in the simplest case, when 
Rashba coefficient is zero. In this case, the diffuson ma- 
trix element I zz takes the form 



^(l^T)2+^ o r2(l + 2(f^-2f) 
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Normalized damped decay rate 




2a k f x 

FIG. 1. The uniform spin dynamics from the weak to the 
strong spin-orbit coupling regime in the presence of both 
Rashba and linear Dresselhaus terms, (a) The normalized 
exponential decay rate Q(Qt) is shown as a function of nor- 
malized Rashba and linear Dresselhaus SOI. (b) The nonzero 
normalized oscillatory frequency, t, , is nonzero whenever 
2ak f r > | or 2/3 i fc / r > §. 

where n so = 2(3 1 k f and 5 = 2&(1 - The corre- 

sponding spin decay rate is 

iflr = 1 - 

2 

(13) 

Eq (1 1 1 II 3|) show that the cubic term will increase the 
exponential decay rate and decrease the oscillatory fre- 
quency. To show the effect of the cubic Dresselhaus term, 
the real(imagine) value of the damped oscillatory fre- 
quency when /?3 7^ is divided by the value when /J3 = 0. 
This ratio is ploted in Fig [5] with respect P3/P1 and 2/3\T. 
When |^ < 0.2, the effect of the cubic term is very small 
and can be neglected. In this case, the damped decay 
rate is always equal to w- as long as fi so > \ and the 
oscillatory frequency Q approach tt so when f2 so r ^> 1. 
This provides a reliable way to measure the momentum 
scattering time r. Further, the strength of the linear 
Dresselhaus SOI can be obtained from Eq. [TT] once we 
know r and the oscillatory frequency from the measure- 
ments. These will be discussed in a later section. 

Now, let us choose a = Pi which is a more unique 
case and gives us the persistent spin helix for special q 
valuesi 14 i 15 ' 17 For the uniform spin polarization, the de- 
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FIG. 2. The uniform spin dynamics from the weak to the 
strong spin-orbit coupling regime in the presence of linear f3i 
and cubic /?3 SOI. (a) The normalized exponential decay rate, 
5ft(zf2r) is constant when P3 is zero and slightly larger than | 
when /?3 is nonzero, (b) The nonzero normalized oscillatory 
frequency, ^s(iflr), appear when Q so t > |. 



cay rate of the spin satisfies 

ifir = 1 - 7T2(!Uf, (14) 

where ft so = 2^/a 2 + fi\ kf. The damped oscillation 
mode will happen when f2 so r = 2y/2akfT > y/2/2, say 
2akfT > 1/2 which is the same as the pure Rashba or 
Dresselhauss case. The oscillating frequency in the clean 
limit is V2Q so = 4akf which is the two fold of the fre- 
quency for the pure Rashba or Dresselhauss case. On 
the other hand, as the real part of zf2r is equal to 1 when 
damped oscillation mode appear, the damped decay rate 
is also the two fold of the case of the pure Rashba or 
Dresselhauss. 



IV. SPIN DYNAMICS AND RABI 
OSCILLATION 

Before we discuss the spin dynamics for the nonuniform 
spin polarization system, let us give a physical explana- 
tion of the result we have obtained. We can construct 
a simple physics picture to describe the spin polarized 
wave theoretically. Taking the Rashba SOI for exam- 
ple, we define the eigenstates |<^) to denote the majority 
band and the \cj) b k ) to denote the minority band. The 
spin of the eigenstate of the SOC 2DEG lies in the x — y 
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plane. The majority electron has opposite spin to the mi- 
nority electron when they have the same wave vector k. 
As a result, The spin polarization along the z direction 



system corresponding to energy differences Ai and A2 
respectively in Fig. [3J 



\ A = 2ak 






\ aY 2) = a„(i + ^) 
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FIG. 3. The dispersion relation due to the linear Dresselhaus 
SOI. The SOI induces the energy gap Ao = 2/3ifc which is the 
spin precession frequency for the single electron spin. How- 
ever, when the system is excited to be a spin polarization wave 
with wave vector q, the spin polarization along the z direction 
is constructed by the superposition of the two electron with 
wave vectors k and k + q. In this case, the spin precession 
frequency will be Axp) — Ao(l ± where Q = 2m/3i. 

can be obtained by the superposition of the majority and 
minority bands as 



vt, 9 = e 



( e _ £/ ) 2 /4 C r 2 ,, 



V2 
1 

71 



4 [E^ £ ' )2/V -^(i4) + i^})]- (15) 



where A is the normalization coefficient, V>t,g i s f ne wave 
function of the system with positive spin polarization 
along z direction with wave vector q and the function 
e ( e - £ / ) / 4cr restrict the spin polarization electrons only 
in the narrow range ^«ej around the Fermi energy e / . 
The expectation value (i/)f !q \cr z cosq'xltp^^) is nonzero 
only when q' = q which confirms that ip^, q can describe 
the spin polarized wave. The energy difference of these 
two electrons in the first (second) term on the right hand 
side of Eq [fT5"|)is ^1(2) as shown in Fig. ©. There- 
fore, can be treated as a collective two level system 

with two Rabi frequencies f2i(2) = A i (2) . For the uni- 
form spin polarization means q = and there is only one 
Rabi frequency flo = 4p Fig. [3] When the system is very 
clean, our results, Eq. (|llll4j) . show that the spin evolu- 
tion is damped oscillation and the oscillatory frequency 
is the Rabi frequency. It is a little surprising that when 
a — (3, although the SOC gap Ao is not a constant, the 
oscillatory frequency is corresponding to the maximum 
splitting energy 4ak f instead of the average splitting en- 
ergy 2y/2akf. In the weak SOC regime, the disorder 
is so strong that the splitting energy due to the SOI is 
completely submerged in the broadening of the band -. 
Therefore, the spin polarization just exponential decays. 
For the non-uniform spin polarization case, since there 
are two Rabi oscillation frequencies f2i and ^2, we ex- 
pect to have two damped oscillatory modes in the clean 



V. NON-UNIFORM SPIN POLARIZATION 

In the case of the non-uniform spin polarized 2DEG, 
the initial state is a spin wave with wave vector q, the 
momentum k is coupled to k + q which makes the center 
of the Fermi sea be shifted to near q. the average mag- 
netic field is nonzero and the off diagonal elements of 
the diffusion matrix appear to couple the different spin 
component. When only considering the Rashba or lin- 
ear Dresselhaus SOI, our numerical calculation does have 
two kinds of spin dynamical modes which are shown in 

Fig.mnai. 
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FIG. 4. The fast oscillatory mode of the nonuniform spin 
dynamics in the strong SOC regime when the system only has 
bulk inversion asymmetry, (a) The normalized exponential 
decay rate, SR(ifir) increase with increasing q and approach 
to one at large q. (b) The nonzero normalized oscillatory 
frequency, 3(if27"), increases linearly at large q, the slope is 
close to 0, 3O t and its value approaches f2 ao (l + jj). 

The damped oscillatory modes at large q can be ap- 
proximately written as 



iflr = 1 



l-(Q so r) 2 (l±|) 2 (16) 



where Q = 2m\i. 

The two damped oscillatory mode and their oscillatory 
frequency satisfy our expectation based on the Rabi oscil- 
lation viewpoint. When q increases, the Rabi frequency 
of the faster mode always increases which makes the 
damped oscillatory mode appear even when f2 so r < i. 
This means we can expect to observe the oscillation for 
the nonzero spin polarization at higher temperature than 
for the uniform spin polarization. In Ref. [ill where 
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FIG. 5. The slow oscillatory mode of the nonuniform spin 
dynamics in the strong SOC regime when the system only 
has bulk inversion asymmetry, (a) The normalized exponen- 
tial decay rate, 5R(ir2r) has a minimum around q — Q and 
approach to one at large q. (b) The nonzero normalized os- 
cillatory frequency, Sj(if2r) is always zero when q is around Q 
and increases linearly at large q. The slope is close to Q 30 t 
and the value approaches fi so (l — q) at large q. 



the spin polarization is uniform, the damped oscilla- 
tory mode appears below 5 K. On the other hand in 
Ref. [T3L where the spin polarization is nonuniform, the 
damped oscillatory mode appears below 50 K. The ma- 
terial, Fermi energy and mobility in these two papers are 
similar. This seems support our Rabi oscillation view- 
point. For the slow oscillatory mode, when q is around Q, 
the corresponding Rabi frequency O2 is around which 
means the spin precession is very slow. Because the Rabi 
frequencies is much smaller than - , the spin polarization 
just decays exponentially and the exponential decay rate 
has its minimum in this regime when q is around Q . A 
particular case is when a = /3i , q = ±4ma and along the 
x direction and (3 3 — defined in Eq.([T]), the Rabi fre- 
quency of the slower mode is zero for all of the electron 
momentum k and Eq. (|16j) becomes the exact solution. 17 
On the other hand, the spin y is a good quantum num- 
ber for all the electron states which means the spin in- 
dependent disorder will never couple the two electrons in 
different bands with different spin directions. Therefore, 
the Rabi frequency of the slower mode is still exactly 
zero even in the presence of the spin independent dis- 
order no matter how strong it is. As a result, the spin 
along the z-direction will never precess and has infinite 
long life time. This provides another way to understand 
the persistent spin helix . 14 i 15 However, the cubic Dres- 
selhaus SOI induces band transition in the presence of 
spin independent impurities and makes the spin life time 
finite. 7 



VI. PROPOSED EXPERIMENTS 

The spin dynamics in the strong SOC regime have sev- 
eral special characters which can be used in experimental 
measurements. 

Momentum scattering time r*: In the spin dynam- 
ics, the Coulomb interaction plays an important role 
in determining the momentum scattering time t * 22 ' 23 . 
This is quite different to the charge transport case where 
electron-electron(e-e) interaction will not change the en- 
semble momentum scattering t p which determines the 
electron mobility. This difference is called spin Coulomb 
drag (SCD). In previous experimental work, SCD was ob- 
served through the spin diffusion coefficient D s — ^VjT* 
by fitting the spin decay rate in the weak SOC regime. 
Here, we provides a way to observe SCD in the strong 
SOC regime by directly measuring the momentum scat- 
tering time t*. Based on Eg 111113"! when only Dressel- 
haus SOI is presented, the damped decay rate is always 
almost equal to 1/2 as long as |^ < 0.2 which is easily 
realized in experiments! 12 > 15 

The strength of SOIs: Here, we would like to empha- 
size that 2(3ikfT — ^ is a very important case and is 
corresponding to the transition point between pure expo- 
nential decay mode and damped oscillatory mode. The 
decay rate at this point is not only equal to but also 
equal to 2 plk f wnen = 0. This means that at this point 
we can obtain the strength of linear Dresselhaus SOI from 
the spin polarization decay rate. When 2/?ifc/r = s, we 
can increase the Rashba SOI by adding a gate voltage. 
As long as < a < pi, according to Eq. [TJJl the spin 
evolution is still decay exponentially and the decay rate 
is (1 - ^Vl -2(fi so r) 2 )/r where n so = 2^u 2 + ffh f 
which gives us the strength of Rashba SOI. 



VII. CONCLUSION 

We have discussed the spin dynamics in the strong 
spin-orbit coupling regime. We describe quantitatively 
the special characters of the damped oscillatory mode 
in this regime. We also compare our result to the pre- 
vious experimental data and find they match very well. 
Based on our theoretical results, a reliable way is pro- 
posed to measure the Rashba and Dresselhaus coefficients 
and electron momentum scattering time which is not cor- 
responding to the mobility due to the Coulomb interac- 
tion. Furthermore, we find that the spin dynamics in 
the 2DEG can be treated as a collective two level sys- 
tem. This helps us semi-quantitatively understand the 
spin dynamics in the strong spin-orbit coupling regime. 
For the nonzero spin polarization case, we predict that 
there exist double damped oscillatory modes at large q 
and explain the persistent spin helix mode from the Rabi 
oscillation point of view. 

We acknowledge Chia-Ren Hu, Ar. Abanov, Yang Liu 
and Victor Galitski for very helpful discussion and sup- 
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Appendix A: Spin dynamic matrix for the uniform spin polarization 

In this section, we derive the spin evolution mode of the uniform spin polarization. According to Eq. [5] the strength 
of SOI is angle dependent and can be written as 

h so = J a 2 +/3ffe/Jl + cos2^cos2(9 + (2(^) 2 - ^ sin(</> + n/4))(l + cos40) - cos(^ + tt/4) cos 2(9 (Al) 
» V A A A 



where costp — Xi/y/ X\ + Af. 

First we consider the case for /?3 = 0. The Hamiltonian is written as 

k 2 

H = — + (a + P)k x a y - (a - f3)k y cr x , 
k 2 

= h Xik x a y + \2k y a x , (A2) 

where k x is along the [110] direction, Ai = a + f3 and A2 = —(a — 0). The Green's function for this Hamiltonian takes 
the form 

E -&i ± 17 + ( a + P) k x°v -(a- f)k y o x 



qR(A) — 2m -*~ 2t 

" ( E ~ L ± i-r) 2 (« 2 + + Jtf cos2 °) 

F, - Jl + J- 4- A, fc-n-.. 4- Aok.rr.„ 

(A3) 



( E - L ± ^) 2 - IM r^ 2 (l + 7cos20) ' 



where r is the momentum scattering time, 7 = 2o fm — T^rri- It is more convenient to write down the element of 



the 2x2 Green's function Eq. (|A3p as 

Gil = G22 = — ( P> ■ r -\- 77; — ^ r 

2 S-|i-AfcVl + 7cos2^±^ ^-H + AWl + 7Cos2^±i 



1 1 1 \/2(— i cost/ 1 cos + sin ^ sin ( 



2 E ~lk~ AfcVl + 7Cos2^± i £ - |i + Afc^/I + 7 cos 20 ± \' VI +7 cos 20 

1 . 1 1 . V2(* cos ^ cos + sin ip sin 

21 " 2 E-^ - AfcVT+7Cos26»± I £ - || + Afc^ + 7 cos 20 ± \ VI + 7 cos 20 

where A = + ^l)/ 2 = V® 2 + P 2 i cos ip = \i/^JX{ + X 2 , and 7 = cos2?/>. 



According to Eq.®, the diagonal element of the spin polarization along z direction has the form 

I zz = ^11,11 ~ ^11,22 ~ ^22,11 + ^22,22 

1 r d 2 k 

~ 2m~^ J (2tt) 2 

The first term and the fourth term in Eq. IA8I are equal to each other and have the form 

1 r d 2 k A R 1 r d 2 k 1 1 1 



fj^_ (r A r R r A r R r A r Ft,, r A r R\ , a r \ 



2mr7 (2tt) 2 11 11 2m J (2tt) 2 4 v £ - e+(fc) - £ 

1 1 

( £ + ft-e_(fc) + 27 + £ + + 

I f^d^ fc+ fc" k+ 



Wmn J v f 1-ittr 1 - iQr + 2i\k^/l + 7 cos 20 1 - iQt - 2i\k^/l + 7 cos 20 1-iQr' 
1 /" 2,r .„, 1 1 1 1 



16tt X d6 ' ( 'l - ifir + 1 - iflr + 2i\k VI + 7 cos 20 + 1 - iOr - 2iAfe VI + 7 cos 26* + 1 - fflr' ) ' 



(A6) 
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where vj = ^f- In the polar coordinate, J d 2 k = J d(k 2 /2)d0. As we assume that \kj <C Ef, d(k 2 /2) ~ mdE where 
m is the effective mass. 

The other two terms are also equal to each other and can be written as 



1 



d 2 k 



2mr J (2?r) 



Lt 21 Lt 12 



1 



1 



d 2 k 1 

j2^4 { E-e + (k)-^ ~ E-e-(k)-£-- 
1 1 



E + {l-e + (k) + £- E + n-e_{k) + £ B 



de 



k+ 



+ 



k- 



16to7t Jq Vf 1 — iQ,T 1 — iQr + 2i\ky/l + 7 cos 20 1 — iftr — 2i\ky/l + 7 cos 20 1 — iClr 



d0{ 



1 

167r 



Substituting Eq. \M\ JA7| to Eq. £H1 we have 

271 



1 - ifiT 1 - zftr + 2zAfcVl + 7 cos 26* 1 - ifir - 2iAfcVl + 7 cos 26> 1 - ifir ' 



(A7) 



1 

27 



5/, d9( 

2tt 

d0 



1 



1 



1 - iflr + 2i\k^/T + 7 cos 20 1 — ifir — 2i\k\/l + 7 cos 20 ; 
1 - iftr 



(1 - iQr) 2 + (0 5O r) 2 (l + 7 cos26l) 



1 - ifir 



dx- 



2tt((1 - zftr) 2 + (0,. t) 2 ) 7 (1 + a cos(x)) ' 



where x = 20 and a = 7(O so r) 2 /((l — i£lr) 2 + (il so r) 2 ). The indefinite integral / dx 1 , * 
Therefore we have 



(A8) 



2arc tanh 



-(-a cos(:c) 



V-l+a 2 



J* 



1 - ifir 



dx- 



2tt((1 - iflr) 2 + {Q so t) 2 ) J (1 + acos(x)) 

(—1+a) tan 

//in- 1 . ; 1 1 1 1 1 

1-iQ.T 



2tt((1 - lilr) 2 + (n so T) 2 ) 
1 - iilr 



2(- 



2arctanh 



(-1+a) tan[§] 
V-l+a a 



V-l + a 2 



V-l + a 2 



2vri 



1 - iflr 



2tt((1 - zftr) 2 + (tt so r) 2 ) V-l + a 2 ^((1 - zftr) 2 + (fi so r) 2 ) 2 - 7 2 (^ t) 4 ' 
When Rashba SOI is zero, the strength of SOIs takes the form 



Plk f Jl + (2^| - 2^)(1 + cos 40). 



Pi P 



(A9) 



(A10) 



To obtain the spin diffusive matrix element J 22 , it is easy to prove that we only need to replace the term Xk^l + 7 cos 20 
in Eq. IA9I with h so in Eq. IA10I Therefore, we have 



r 



2(1 - ifir) 



4tt / ( (1-^t) 2 + (2/i so t) 2 ' 

1 - ifir 



r)d0 



^(l- l nr) 2 + n 2 so r 2 ^(l - ^^7r) 2 + fi 2 r 2 (l + 2(f) 2 - 2f ) 

I 
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